Design of reflective, photonic shields for atmospheric reentry 
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We present the design of one- dimensional photonic crystal structures, which can be 
used as omnidirectional reflecting shields against radiative heating of space vehicles 
entering the Earth’s atmosphere. This radiation is approximated by two broad bands 
centered at visible and near-infrared energies. We applied two approaches to find 
structures with the best omnidirectional reflecting performance. The first approach 
is based on a band gap analysis and leads to structures composed of stacked Bragg 
mirrors. In the second approach, we optimize the structure using an evolutionary 
strategy. The suggested structures are compared with a simple design of two stacked 
Bragg mirrors. Choice of the constituent materials for the layers as well as the 
influence of interlayer diffusion at high temperatures are discussed. 
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I. INTRODUCTION 


Reflecting photonic systems (RPS) are of high interest for many practical applications 1 , 
including potentially as protective coatings for atmospheric reentry. During atmospheric 
reentry, space vehicles experience extreme environments which put high performance de- 
mands on their thermal protection systems. A major source of heating during reentry is 
radiation generated in the shock layer. Extreme temperatures ionize gaseous species, pro- 
ducing radiation which impinges on the surface of the vehicle. This radiative heating scales 
like v 8 and can be a major factor especially in manned mission from the Moon or Mars. 
For Jupiter entries, 99% of the heating is radiative. Because planetary atmospheres have 
different chemical compositions, their radiation profiles will also be different. Therefore, an 
ideal RPS solution will reflect radiation in specified wavelength ranges for a given planet. 
The schematic illustration of the thermal protection system (TPS) of a vehicle is shown in 
Fig. 1. 
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Figure 1: Schematic of a possible thermal protection system for a space vehicle. RPS protects 
against radiative heating. Thermal barrier coating retard conductive heat propagation. 

Shock layer radiation spectra at peak heating were modeled and analyzed 2-4 . Experimen- 
tal data on radiation spectra of shock- heated air were obtained at flight-relevant conditions 5 . 
Radiation spectra were measured in the range of 190-1210 nm. These spectra contain two 
broad peaks which are mostly due to radiating species of N 2 , Nj, N and O. The first peak is 
approximately in the range of 650-1200 nm (B{) and the second is in the 300-500 nm range 

m. 

The goal of this paper is to propose an optimal design of a one-dimensional photonic 
structure with high omnidirectional reflection in these two wavelength intervals. The term 
“optimal” implies the best reflecting performance of the structure for a fixed number of 
layers. This design principle is quite general and can be used to make omnidirectional 
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RPS for other wavelength ranges as well. The advantage of the proposed one-dimensional 
structures is that they are relatively simple to model and to fabricate. 

The constituent materials silicon dioxide (SiC^) and zirconium dioxide (Zt0 2 ) are cho- 
sen. These materials meet a number of the strict requirements of this application, such as 
high-temperature stability, low absorptivity, oxidation resistance and others. However, the 
Si02/Zr02 couple has small dielectric contrast, making the design of a omnidirectional RPS 
more difficult than for materials with a higher dielectric contrast. 

This paper is organized as follows. In Section II, the optimization model of an omnidi- 
rectional RPS appropriate for two energy intervals is defined. Next, in Section III, we apply 
a band gap analysis to Si 02 /Zr 02 periodic multilayers and construct an omnidirectional 
RPS in the desired energy intervals. By combining (stacking) different periodic segments in 
one structure a high reflectivity can be achieved. In Section IV we briefly discuss different 
optimizers, which are applicable for the optimization of layered structures. As an example, 
we optimize the omnidirectional reflectivity of a layered structure using ES and compare our 
results with the Bragg and optimized Bragg mirrors. Then, in Section V optimization of the 
RPS for either to B\ or to B 2 energy interval is performed. The optimization of the RPS for 
both energy intervals is performed in Section VI. Finally, in Section VII changes of in the 
multilayer structure due to high temperature interdiffusion and the possible consequences 
on reflecting properties of RPS are discussed. 


II. DEFINITION OF THE PROBLEM 

The RPS for Earth reentry should be optimized for energy intervals Bi and B 2 and 
should possess omnidirectional reflectivity. The energy ranges for B\ and B 2 are 1-2 eV 
(~ 650 — 1200 nm) and 2.4-4. 2 eV (~ 300 — 500 nm) respectively. We have chosen sili- 
con dioxide (Si 02 ) and zirconium dioxide (Zr 02 ) as the materials for the individual layers. 
Silicon dioxide and zirconium dioxide have high melting temperatures T m = 1650 °C and 
T m = 2715 °C respectively and are optically transparent in the ranges of interest. Si0 2 
(quartz) and ZrC >2 have similar coefficients of thermal expansions. The high-temperature 
interaction of SiC >2 and ZrC >2 layers results in the formation of a high-temperature stable 
phase of ZrSiCR which may have possible diffusion barrier properties (see Section VII) . This 
phase may prevent further Si and Zr interdiffusion, and therefore result in reduction of pho- 
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Table I: Melting points (T m ), dielectric permittivity (i), and density (p) of some possible 

materials for the high temperature photonics 


Material 

T m , °C 

e, Ref. 6 

p, g- cm 3 

MgO 

2825 

3.7 + 0.6i (at 1 pm) 

3.6 

y 2 0 3 

2440 

3.6 (at 1 pm) 

5.0 

A 1203 

2050 

2.7 (at 1 pm) 

4 

Ti0 2 

1840 

7.5 (at 1 pm) 

4.2 

Si0 2 

1650 

2.1 (at 1 pm) 

2.2 

Zr0 2 

2715 

4.8 (at 0.8 pm) 

5.9 

Hf0 2 

2800 

3.6 (at 0.8 pm) 

9.7 

MgF 2 

1260 

1.9 (at 0.8 pm) 

3.1 

SiC 

2830 

6.8 (at 1 pm) 

3.2 

Si 3 N 4 

1900 

4 (at 1 pm) 

3.2 

Nb 2 0 5 

1500 

4.4 (at 1 pm) Ref. 7 

4.5 


tonic properties. Permittivity data of the materials are shown in Fig. 2. Both materials 
are dielectrics with a wide gap, so we assume that permittivity should not strongly change 
with temperature. However, all calculations are performed with permittivity values at room 
temperature, since data for permittivity at high temperature are not available. Basic prop- 
erties for possible candidate materials for high-temperature photonic applications are given 
in Table I. 
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Figure 2: Permittivity of SiC >2 and ZrC >2 from 6 . 


Our goal is to design a layered structure which maximizes the percentage of reflection 
R of unpolarized electromagnetic radiation within two energy intervals for a fixed number 
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of layers N. Therefore, we define a fitness function as the average reflection in two energy 
intervals: 


Fitness = {R(di,d 2 

<*>;<*> <i, 

(1) 

1 1 m t 

cj, 9) + Rp(ui, 9))dud9 



f f dudO 

(2) 


^ 1,2 0 


where R s and R p are the reflection coefficients for s- (electric field parallel to layers) and 
p-polarized (magnetic field parallel to layers) light respectively and are integrated over the 
angles of incidence 9 E [0, |] and energy intervals hui E [1, 2], hcu 2 E [2.4, 4.2] eV. The 
thickness of the i-th layer (i = 1, ..., N) is d*. 

We need to maximize the function R(d\, d 2 , ... , dw) of N variables, which are layer thick- 
nesses. Moreover, the obtained solution should also be robust. This implies that reflection 
(R) does not vary strongly with small deviation of layer thicknesses Sdi. In other words, 
unavoidable fabrication inaccuracies of layer thicknesses should not significantly change re- 
flecting performance of the structure. 

III. DESIGN OF THE RPS USING PHOTONIC BAND GAP APPROACH 

Construction of omnidirectional mirrors from layered media were analyzed and discussed 
in many works 8-10 . There are two necessary requirements to achieve omnidirectional re- 
flection. The first is a sufficiently high dielectric contrast for the constituent materials. 
The second is that the dielectric constant of the lower index material is larger than the 
dielectric constant of the ambient medium by some critical amount. Reflection in a broader 
energy interval can be achieved by stacking several periodic layered media with different 
periods 11 . Another way to achieve reflection in a broad range is to use a chirped periodic 
layered medium 11,12 . Other optimization procedures, like the quasi-Newton algorithm 13 or 
genetic algorithms 14,15 have also been applied to improve the performance of finite multilayer 
mirrors. 

Let us consider an inhnitely-periodic structure consisting of Si02 and Zr02 layers. The 
geometry is shown in Fig. 3. a. For normal incidence light {k x = 0) the gap-midgap ratio 
is maximized when d±n\ = d 2 n 2 (quarter wave stack condition). For grazing incidence 
(k x — > cu/c), the gap-midgap ratio is maximized when d\\J n\ — 1 = d 2 \J n \ — l 16 . 
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(a) (b) 

Figure 3: a) Multilayer SiCH/ZrC^ with period A constant through the whole stack, dg 2 and 
are thicknesses and refractive indexes of the layers respectively, b) Multilayer stack consisting of 
segments with different periods. Angle of light incidence is 9, k = (k x ,0,k z ) is the wave vector. 

The main difficulty in designing a RPS with sufficiently broad band gaps for both Bi 
and B 2 energy intervals is the low dielectric contrast of SiC >2 and ZrC >2 materials. A simple 
periodic Bragg mirror is not a suitable solution, because the band gaps are narrow and not 
omnidirectional. The band structure of a Bragg mirror with period A = 425 nm is shown 
in Fig. 4. The band gap for p-polarized light becomes narrower for the inclined incidence 
light as can be seen in the plot. Since point B is slightly above point C, the structure has 
110 omnidirectional reflection. 

In order to circumvent the difficulty of narrow gaps, we can stack segments with different 
periods, as shown in Fig. 3.b. If the periods of the segments are chosen appropriately, 
then the total omnidirectional band gap can include the entire desired energy range. As 
can be concluded from Fig. 4, the condition for band gap overlapping in p-polarization is 
that point B of the next segment with a smaller period should be slightly below point A 
of the previous segment with a larger period. Band gap overlapping in s-polarization is 
automatically fulfilled using this criteria. 

In order to generate an omnidirectional gap in the first range (1-2 eV) we need ten 
periodic segments with uniformly spaced periods from 200 to 425 nm. The band structure 
of ten segments is shown in one of the plots in Fig. 5. a. Note that the B 2 range 2. 4-4. 2 meV 
is already significantly covered with band gaps of higher orders. This means that reflection 
of finite structure for the B 2 range will also be high. However, in order to completely cover 
the B 2 range with first order gaps, we need eight additional segments with uniformly spaced 
periods from 90 to 180 nm. The resulting band structure is shown in Fig. 5.b. Thus, a total 
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of 18 segments are required to cover both B\ and B 2 energy intervals with band gaps. 

The magnitude of Im(fc z ) in Fig. 5 determines the rate of exponential decay of the prop- 
agating wave and also affects the reflectivity of the structure. As seen in the plot, the 
magnitude of Im(fc z ) is larger for higher energies (or smaller periods of segments) in both B\ 
and B 2 ranges. However, the product Im(fc z )A remains constant for all periods A in the seg- 
ments. Therefore, we take equal numbers of layers in each of the 18 segments. Calculations 
of the reflectivity of these stacked segments are discussed in Sections V and VI. 



Figure 4: Band structure of Si02/Zr02 Bragg mirror with period 425 nm for s- and p-polarized 
light. White region between shadowed regions corresponds to band gap. Solid black line denotes 
air light cone, dashed and dash-dotted lines denote Brewster’s cones of Si02/Zr02 and air/Zr02 
interfaces respectively (ZrC >2 is on top of structure). 


IV. ES OPTIMIZER FOR LAYERED STRUCTURES 

Until now, our approach to designing structures with good reflection was based on band 
gap analysis. We stacked 18 segments of Si0 2 /Zr0 2 with different periods in one structure 
to obtain reasonable performance. This is one possible way to design. However, it is not 
necessarily the best one and it is interesting to find out how close we are to the optimum. 
Different optimizers 17,18 can be helpful to find good solutions, that are often difficult to find 
analytically. 

If the structure consists of N layers, then there are N parameters (thickness of each 
layer) for the optimization. The search space should be defined for each parameter within 
a reasonable interval [rf min , d max ]. For such optimizations different evolutionary algorithms, 
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Figure 5: a) Band structure of ten segments with uniformly spaced periods from 200 to 425 nm 
shown in one plot. The total omnidirectional band gap is in the B\ energy interval 1-2 eV. b) 
Band structure of eight segments with uniformly spaced periods from 90 to 180 nm shown in one 
plot. The total omnidirectional band gap is in the B 2 energy interval 2.4-4. 2 eV. Colorbar 
corresponds to imaginary part of k z in pm -1 units. 


like genetic algorithm (GA), micro-genetic algorithms (MGA), or evolutionary strategy (ES) 
may be applied. The optimization task becomes difficult, if there are many parameters for 
optimization, and/or if the fitness function has a complicated behavior with many local 
optima. We used ES for the optimization of our RPS. In previous applications ES usually 
outperformed GA and MGA. 

As an example of the ES approach, we maximized the average reflection ( R ) of the 
Si02/Zr02 structure containing N=8 layers, at wavelength A = 0.8 pm. Permittivities of 
Si02 and Zr02 are taken to be Ei = 2.1 and £2 = 4.8 at A = 0.8 pm respectively. Layers are 
embedded between air interfaces. 

The average reflection to be optimized is defined as: 

(R) = R(d 1 ,d 2t ... ,d„)= f ’ /2 Ml±M!l d o, ( 3 ) 

Jo 7r 

where R s and R p are reflectivities for s- and p-polarized light respectively, di is the individual 
layer thickness, and 6 is the angle of light incidence. The results of ES optimization are shown 
in Fig. 6. a. Dash-dotted curve corresponds to a Bragg mirror with layer thicknesses: 

, A , A 

“1 — , . — j “2 — , , — • 

4-v/^l 4y ^ 
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( 4 ) 
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Figure 6: a) Dash-dotted curve - reflectivity of the Bragg mirror, dashed curve - reflectivity of 
the optimized Bragg mirror, solid curve - reflectivity of the optimal structure. Wavelength 
A = 0.8 /mi b) Thicknesses of individual layers for three structures, number of layers N=8. 


Dashed curve in the plot corresponds to an optimized Bragg mirror, with two optimized layer 
thicknesses. Solid line corresponds to a structure where all eight thicknesses are optimized. 
As seen from the plot, the optimization of all layer thicknesses (solid line) leads to the 
best performance. This reflects the fact that a Bragg mirror is not an optimal structure 
for omnidirectional reflection. A comparison of the layer thicknesses for three structures is 
shown in Fig. 6.b. 


V. OPTIMIZATION OF THE RPS ADJUSTED TO ONE ENERGY 
INTERVAL 

First we performed optimization of the RPS with respect to only one energy interval. 
Let us first consider RPS for the B\ energy interval. The fitness function to be optimized 
is defined in the expr. (2). Band gap analysis revealed that 10 stacked periodic segments 
may cover this range with band gaps. The dashed line in Fig. 8. a shows reflectivity ( R\ ) as 
a function of the number of layers in this structure. In order to improve the reflectivity, we 
performed an ES structure optimization in two ways. 

In the first way, optimization was done in the full search space. This means that the 
thickness of each layer could have any value below d max = A max /2n, A max = 1.2 /im. This 
yields d max = 280 nm (ZrC^) and d m ax = 425 nm (SiCL). As mentioned earlier, optimization 
becomes demanding, if there are too many parameters for optimization and/or if the fitness 


9 


function has a complicated behavior with many local optima. Therefore, it is helpful to 
visualize the dependence of the fitness function on the optimized parameters. Fig. 7 shows 
the dependence of (Ri) (fitness) on the thicknesses of two arbitrary layers in the structure (4- 
th and 10-th layers are taken) containing 20 layers. One can see that the fitness function has 
a rather smooth behavior, which is advantageous for optimization. Moreover, smoothness of 
the fitness function means that the obtained solutions are quiet robust and small deviations 
(fabrication inaccuracies) in layer thicknesses will have a small impact on the reflectivity. 
Results of the optimization in the full search space are shown in Fig. 8. a (dash-dotted line). 



' d 10 , nm 

Figure 7: (Ri) (fitness) of the structure as function of the thicknesses of two layers (d^ and dio)- 

Structure consists of 20 layers of Si02/Zr02- 

In our second optimization approach, we reduced the search space by using the thicknesses 
( d s i e9 ) of 10 stacked periodic segments as a first order approximation for the ES optimizer. 
Then we varied the thickness of each layer in the range of Sd = [—50, 50] nm to obtain the 
optimum fitness. Thus, the obtained optimal thicknesses can be represented as: 

d° pt = df 9 + Sdi, i = 1, . . . , N. (5) 

Results for the optimization in the reduced search space are shown in Fig. 8. a (solid line). 
In both optimization procedures we restricted the number of fitness calls to 120000 in order 
not to make the optimization very time-consuming. Of course, the reflectivity of the RPS 
can be further improved by means of simply increasing the time of optimization. Therefore, 
it should not be confusing that for N=40 in Fig. 8. a dash-dotted line is below the solid line. 

As seen in Fig. 8. a, the structure optimization in the full search space (dash-dotted curve) 
becomes less efficient (time-consuming) for a high number of parameters N > 20 and the best 
performance comes from the optimized stacked segments. Thicknesses of individual layers 


10 



Figure 8: a) Reflectivity (Ri) as function of the number of layers N in the structure. Dashed line 
- stacked periodic segments (periods from 200 to 425 nm), solid line - stacked segments optimized 
with ES, dash-dotted line - solution of the ES optimizer in the full search space, b) Thicknesses 
of individual layers of the three structures, total number of layers is 20. 

for the three structures are compared in Fig. 8.b. The improvement in reflectivity, however, 
is rather small and does not exceed 3%. We can therefore conclude that the stacked periodic 
segments is already a good solution for high reflection in B\ and only minor improvements 
can be achieved. 

A similar procedure was also performed for the second energy range B 2 . Results are 
shown in Fig. 9.b. 

VI. OPTIMIZATION OF THE RPS ADJUSTED TO TWO ENERGY 
INTERVALS 

In a similar way we can optimize the reflection of RPS with respect to two energy intervals 
simultaneously. Now the fitness function to be optimized is ( R ) in the expr. (1). Band gap 
analysis revealed that 18 stacked periodic segments may cover B\ and B 2 intervals with 
band gaps. Again we can apply two approaches to perform the structure optimization using 
ES. The first approach is an optimization in the full search space, where the thickness of 
each layer is any value below d max = A max /2n, A max = 1.2 /im. However, as we saw in 
Section V, this optimization procedure is not very efficient for many parameters N > 20. 
So we used only the second approach, which optimizes 18 stacked segments in the reduced 
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Figure 9: a) Reflectivity (R 2 ) as function of the number of layers N in the structure. Dashed line 
- stacked periodic segments (periods from 90 to 180 nm), solid line - stacked segments optimized 
with ES, dash-dotted line - solution of the ES optimizer in the full search space, b) Thicknesses 
of individual layers of the three structures, total number of layers is 16. 



Figure 10: a) Reflectivity (R) of the structure for B 1 and B 2 energy intervals as function of the 
number of layers. Dashed line - stacked periodic segments, solid line - stacked segments optimized 
with ES. b) Thicknesses of individual layers of the two structures, total number of layers N=36. 

search space. We varied the thickness of each layer in the range [-30, 30] nm to And the 
optimum fitness. The obtained results are shown in Fig. 10. As can be seen, optimization 
of 18 stacked segments (solid line) improves the reflecting performance. The improvement, 
however, is not large ~ 3% for 36 layers and ~ 0.5% for 72 layers. So, if the structure 
contains more than 36 layers, stacked periodic segments are close to the optimal solution for 
high reflectivity in the Bi and B 2 energy intervals. 
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For the sake of completeness, it is also worth plotting reflectivity spectra of the obtained 
structures. Up to now we plotted only the average values of reflectivity. Reflectivity as a 
function of energy of the stacked segments (dashed line) and optimized stacked segments 
(solid line) is shown in Fig. 11. At each energy point reflectivity is averaged over angle of 
incidence 9. Dash-dotted line in Fig. 11 is shown for the sake of comparison with the simple 
design, where two Bragg mirrors are stacked in one structure. Band gaps of these Bragg 
mirrors are centered in the middle of B i and B 2 intervals, the corresponding periods are 
Ai = 235 nm and A 2 = 105 nm. A single Bragg mirror will partially cover only one interval 
(B\ or B 2 ). When Bragg mirrors are stacked together, both intervals are partially covered. 
Therefore, two peaks are observed in Fig. 11. It is obvious that optimized segments (solid 
line) show the best reflecting performance. 



Figure 11: Reflectivity (Re) as function of energy. Dashed curve - 18 stacked periodic segments, 
solid curve - 18 stacked segments optimized with ES, dash-dotted curve - reflectivity of the 
stacked Bragg mirrors with periods Ai = 235 nm and A 2 = 105 nm. 


VII. POSSIBLE INSTABILITIES OF THE DESIGNED PHOTONIC 
STRUCTURE DUE TO HIGH-TEMPERATURE INTERACTION OF 
INDIVIDUAL LAYERS 

At high-temperatures, the interaction of Zr0 2 and Si0 2 layers (Fig. 12) may be a source of 
dynamic instability in the position and width of the photonic band gap in designed photonic 
heat-shields. At temperature above ~ 1300 °C 19 , irreversible reactions Eq. (6) between of 
Si0 2 and Zr0 2 , will cause refractory zirconium orthosilicate ZrSi 04 (zircon) interlayers to 
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Figure 12: Scheme illustrating formation of ZrSi 04 interlayers as a result of high-temperature 
interaction Si02 + ZrCR — >• ZrSi 04 under external heat radiation in the wavelength ranges B\ and 
2?2- Thickness of intermediate ZrSi 04 layers after exposure during operation time is much smaller 
than thickness of the SiCR and ZrCR layers d\ and (R (5 <C d\, 5 <C cR) due to high activation 
energy of reaction of ZrSiCR formation. That is why erasing ZrSiCR interlayers do not change 
significantly position and width of the photonic band gaps of designed multilayers. 


form at the interfaces SiCR/ZrCR: 


SiCR + ZrCR — * ZrSiCR (6) 

The upper limit of ZrSi0 4 stability lies between 1669 and 1683 °C 20 which is close to the 
lowest eutectic temperature of 1687 °C. Reaction Eq. (6) remains irreversible until 1550 °C 
at which point, a small fraction of ZrSiCR decays into constituent oxides. The rate of the 
ZrSiCR decomposition reaction increases at 1650 °C 19 . The decay of ZrSi0 4 at 1673 °C was 
established based on an analysis of literature data, experimental results and thermodynamic 
assessments 19 . The lowest eutectic in the system was established at the slightly higher 
temperature of 1687 °C 19 . These results are in good agreement with DTA investigation of 21 
where the zircon decomposition was determined to be at 1674 °C and the eutectic melting 
at 1680 °C. The preliminary literature analysis indicates that diffusivity of Si in SiC> 2 22 ’ 23 
is several orders of magnitude higher than diffusivity of Si in ZrSiCR 24 and Zr in ZrC> 2 25 . 
Diffusivity of Zr in ZrSiCR is assumed to be similar to diffusivity of Hf in ZrSiCR, which is 
higher than diffusivity of Si 24 . Among tetravalent cations, the diffusivity of Si in ZrSiCR is 
approximately the same as that of Ti but about one order of magnitude faster than diffusivity 
of Hf and two orders of magnitude faster than diffusivity of U and Th 24 . A low speed of 
reaction Eq. (6) and a small thickness of formed ZrSiCR interlayer after exposure time of the 
photonic heat-shield of approximately 300 sec are expected (see Fig. 5 in 24 , measurement of 
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temperature response for standard PICA sample at 570 W/cm 2 , 0.42 atm, total heat load 
of 14’250 J • cm 2 , surface temperature ~ 2760 °C). Based on diffusion measurements, the 
activation energy of diffusion of Si in zircon was calculated 22 , which is higher in comparison 
with other silicates. These data also confirm slow rates of reaction 6 at 1200 °C indicated 
by 19 . 

The intermediate ZrSiCB layers consist of atoms of both contacted materials, and therefore, 
their dielectric permittivity has an intermediate value between the dielectric constants of 
ZrC >2 and SiC >2 materials. This affects reflection of light at the boundary. The influence of 
the intermediate layer on the transparency of photonic crystals has been considered in 26 . 
It was found that the effect depends on the relation between the light wavelength and the 
interlayer thickness 5. Indeed, to estimate the effect, we can write the dielectric permittivity 
of the layer as a function of distance from the boundary e = e(z) (here e{z) = £i, if 
z < —S, but e(z) = £ 2 , if z < +5). The effect depends on the relation between the 
thickness of the layer (5) and the light wavelength (A). Namely, the commonly used boundary 
conditions for the light wave fields hold at the abrupt boundary S <C A. This leads to the 
reflection and transmission coefficients used in this paper. On the contrary, at the smooth 
boundary 5 S> A we can use the semiclassical approximation and the wave equation with the 
position dependent wavelength A(z) = 2-Kc/uj^/e^z), which leads to the exponentially small 
reflectivity. Transition between these limiting cases for the similar quantum mechanical 
problem has been studied in 27,28 . Thus, the results presented here hold for 5 <C A, i.e. 
if the diffusion length is much smaller than the thickness of the layers we have deposed. 
Apparently, evolution of the interlayer thickness makes the reflectivity of the multilayer 
coating time- dependent. Nevertheless, the expected thickness of formed ZrSiC >4 interlayer 
should not exceed several lattice constants (25 ~ 10 nm) after heat-shield operation time. 
This is significantly smaller than the wavelength of external heat radiation in both B\ and 
B 2 ranges and, therefore, does not change essentially the results of our modeling. It was 
found 26 that such small thickness of the ZrSiCB interlayer can change the depth, but not the 
position and width of the band gap of designed multilayer photonic structure. 
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VIII. CONCLUSIONS 


Modeling layered Si02/Zr02 structure for potential use as photonic shield for Earth reen- 
try was performed both analytically and numerically for shock layer radiation in near-IR and 
near-UV ranges. The analytical approach indicates that high omnidirectional reflectivity can 
be achieved by means of stacking 18 periodic segments in one structure. Optimization using 
evolutionary strategy resulted in structures with improved reflection performance. These 
structures demonstrate significantly enhanced reflectivity than a simple stack of two Bragg 
mirrors. The ES structure optimization method is quite general and can be used to design 
reflectors for other wavelength ranges without any restrictions on the constituent materi- 
als. High-temperature interactions at the Si02/Zr02 interface resulting in the formation of 
ZrSi0 4 interlayers is expected to be slow. Therefore, the interlayer is not expected to have 
a significant impact on the reflecting performance of the suggested multilayer coatings. 
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